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1.0  INTRODUCTION 

The  following  report  summarizes  the  work  done  on  a  one  year  research 
project  on  microwave  acoustic  devices,  principally  topics  related  to  the 
thin  film  resonator,  TFR. 

Prior  and  concurrent  work  supported  by  AFOSR  was  on  thin  film  materials 
of  ZnO,  AlN,  and  LiNbO^  deposited  by  the  sputtering  process.  These  films 
were  evaluated  in  the  TFR  device  configuration.  In  the  course  of  the 
materials  program,  it  became  apparent  that  the  TFR  is  a  potentially  valuable 
device  per  se  whose  device  physics  problems  need  further  study. 

Device  physics  problems  of  the  TFR  that  required  further  study  are  the 
multidimensional  boundary  value  problem  of  microwave  acoustic  modal 
analysis,  temperature  and  aging  effects,  and  limits  of  TFR  performance  in 
the  1-2  GHz  frequency  range. 

Because  of  the  limited  duration  of  this  project  (one  year),  the  topics 
had  to  be  limited  to  the  personnel  and  prior  experience  resources  available. 
In  particular  the  two  dimensional  numerical  analysis  of  the  TFR  had  been 
formulated  but  needed  further  refinement.  A  prior  project  from  RADC  that 
had  initiated  a  study  of  doping  effects  on  acoustical  properties  in  GaAs  was 
followed  here  by  a  study  of  deep  Zn  diffusions.  Finally,  a  study  of 
acoustical  mode  trapping  effects  was  initiated  in  order  to  determine  the 
lower  limits  of  resonator  Q  in  the  1-2  GHz  frequency  range.  Because  of  the 
limited  duration  of  funding  only  MSEE  degree  level  efforts  could  be 
supported. 

The  next  three  sections  is  a  summary  of  the  projects  and  progress  made 
during  the  one  year  of  funding. 


2.0  HIGH  Q  RESONATORS 


2.1  Introduction 

Resonator  Q  is  an  important  parameter  that  characterizes  resonator 
performance  in  oscillator  and  filter  circuits.  The  resonator  Q  is  related 
to  the  quality  of  the  material  and  topology  of  the  films. 

The  thin  film  resonator  is  a  material  composite  in  the  configurations 
shown  in  Figs  2.1  and  2.2.  Research  efforts  were  directed  toward  improving 
the  Q  of  the  Si  base  resonators  by  adjusting  sputtering  conditions  to 
improve  film  quality  and  electrode  design  in  the  configuration  of  Fig  2.1. 


2.1.1  Film  Substrate  Evaluation 

The  most  definitive  piezoelectric  film  evaluation  is  the  measurement  of 
the  film's  piezoelectric  coupling  coefficient  and  Q. 

Measurements  of  unloaded  0  are  inferred  from  the  resonator  impedance  as 
derived  from  transmission  line  reflection  coefficient  data  measured  on  an  HP 
8510T  network  analyzer.  The  resonator  Q  is  calculated  from  the  measured 
impedance  data  using  the  relation 


Q 


fr  dZ» 

2  df  fr 


where  Z+  is  the  phase  of  the  impedance  and  fr  is  the  resonant  frequency 
defined  implicitly  as  the  frequency  of  maximum  Q. 

The  electromechanical  coupling  coefficient  K2  of  a  resonator  is  also  an 
important  quantity  for  filter  synthesis  applications  because  it  affects  the 
filter  performance  in  two  ways.  First,  the  pole-zero  separation  is 
dependent  upon  K2  with  larger  K2  being  capable  of  synthesizing  wider 
bandwidths.  Second,  the  resonant  series  resistance  is  scaled  by  K2  and  the 
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capacitive  reactance  of  the  structure.  For  a  one  dimensional  resonator 
which  has  a  very  large  diameter  to  thickness  ratio,  the  coupling  coefficient 
is  determined  from 

K2  -  *r/tan*r 

where 

♦r  -  I  fs/fP 

fs  -  series  resonant  frequency 
fp  -  parallel  resonant  frequency 

The  K2  determined  from  this  relation  for  nonideal  resonators  is 
referred  to  as  an  effective  K2  since  the  value  is  not  simply  related  to 
fundamental  material  parameters. 

The  best  method  of  measuring  these  macroscopic  film  properties  is  to 
use  a  resonator  in  one  of  the  accepted  configurations.  In  all  such  cases 
the  measurement  is  subject  to  artifacts  due  to  the  measuring  structure. 

These  macroscopic  factors  are  discussed  below. 

Energy  Trapping.  The  degree  of  energy  trapping  is  a  function  of  the  device 
geometry  and  specific  resonant  mode.  Resonator  Q  is  a  function  of  energy 
trapping  and  measurements  indicate  a  minimum  material  Q  for  a  specific  mode. 
Thus,  a  resonator  structure  that  allows  transverse  acoustic  radiation  along 
the  membrane  and  into  the  substrate  would  exhibit  a  lower  Q. 

Electrode  Metallization  Effect.  Resonator  acoustic  losses  are  increased  by 
the  relatively  high  acoustical  loss  in  metals  compared  to  the  piezoelectric. 
Thinner  metal  electrodes  reduce  mechanical  losses  but  increase  the 
electrical  resistance.  Note  that  increased  metal  thickness  immediately 


within  to  the  resonator  will  cause  the  wave  to  migrate  to  the  thicker  or 
slower  structure  and  hence  lead  to  trapping.  The  electrode  resistance 
effect  is  most  pronounced  at  series  resonance  where  large  displacement 
currents  in  the  piezoelectric  must  flow  through  the  electrode  metal.  This 
affect  is  observed  as  a  lower  series  resonant  Q  as  compared  to  parallel 
resonance  Q. 

Geometrical  Defects.  Surface  contamination  reduces  resonator  Q  and  causes 
frequency  shifts.  Etching  imperfections  such  as  Pits  or  hillocks  give  false 
resonator  operation  due  to  mode  shifts  and  scattering  losses. 

Stray  Capacitance.  Stray  capacitance  in  the  electrodes  appear  in  parallel 
with  resonator  and  causes  the  parallel  resonant  frequency  to  decrease.  The 
reduced  pole-zero  spacing  reduces  the  apparent  coupling  coefficient  and 
lowers  the  achievable  filter  bandwidth. 

Interpretation  of  Q  and  K  measurements  must  be  done  with  some  caution. 
The  K2  as  defined  by  the  phase  slope  definition  of  series  and  parallel 
resonance  may  depart  from  the  material  K2  for  low  Q  resonators. 

Consequently,  a  change  in  material  growth  conditions  or  resonator  geometry 
that  appears  to  have  increased  K*  may  in  reality  have  reduced  Q  and  the 
resultant  spread  in  phase  incorrectly  interpreted  as  increased  material  K  . 

In  order  to  evaluate  films  with  the  least  ambiguity  we  have  adopted  the 
following  procedures: 

1.  Films  not  in  resonator  configurations  or  deposited  over  the  whole 
wafer  are  optically  examined  by  Normarski  microscopy,  x-ray  diffraction,  and 
when  appropriate  Auger  chemical  analysis  (to  determine  oxygen  in  AlN  to  0.1% 
levels)  or  electron  microscopy. 
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2.  Resonator  film  structures  are  first  examined  by  Normarski  phase 
contrast  microscopy  to  find  structural  defects  and  subsequent  electrical 
evaluation  is  done  by  one-port  microwave  reflection  measurements.  From  the 
electrical  data,  series  metallization  resistance,  series  and  parallel 
resonant  frequencies  and  respective  Qs,  and  finally,  the  uncorrected 
piezoelectric  coupling  coefficient  are  obtained. 

2.1.2  Resonator  Analysis  and  Modeling 

Resonator  analysis  and  modeling  takes  two  basic  forms.  These  are  the 
one-dimensional  models  used  for  ideal  resonators  and  the  approximate  lumped 
element  circuit  models  (Butterworth  Van  Dyke)  used  for  some  synthesis 
applications.  Only  the  one  dimensional  Mason  model  (Fig.  2.3)  can  be 
derived  from  first  principles.  A  two  dimensional  numerical  analysis  based 
model  is  described  in  section  3.0. 

The  one  dimensional  Mason  model  is  extremely  useful  because  it  can  be 
used  to  model  the  performance  of  a  single  mode,  and  with  modification, 
multimode  two  and  three  dimensional  resonators.  Basically,  the  Mason  bulk 
wave  model  treats  the  piezoelectric  plate  as  a  three  port  network  having  one 
electrical  port  and  two  acoustical  ports. 

2.2  Experimental  Results 

The  measured  results  shown  in  Fig.  2.4  were  for  an  AlN/p+Si  membrane 
resonator  fabricated  by  selective  etching  of  the  Si  membrane  followed  by  dc 
planar  magnetron  deposition  of  the  AlN.  Series  resonance  Q  was 
approximately  1500  and  parallel  resonance  Q  was  5000.  The  decreased  series 
Q  was  found  to  be  due  to  a  finite  (one  ohm  per  square)  sheet  resistance  in 
the  metal  electrode. 
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However,  energy  trapping  was  still  accomplished  by  the  electrode 
structure  as  evidenced  by  the  high  Q  and  multitude  of  anharmonic  resonances . 
Thus  the  research  goal  of  demonstrating  energy  trapping  was  accomplished. 


3.0  MASKED  DEEP  ZINC  DIFFUSIONS  IN  GaAS 

3.1  Introduction 

Temperature  compensation,  TC,  of  thin  film  resonators  is  an  important 
feature  for  eventual  device  applications.  In  composite  resonators  TC  can  be 
accomplished  by  a  positive-to-negative  balance  between  two  or  more  layers  in 
the  structure.  This  occurs  in  the  composite  resonators  using  p+  silicon  as 
a  support  structure.  The  experimental  results  1  and  theory^  for  this  case 
are  well  documented. 

The  objective  of  this  task  was  to  evaluate  similar  effects  in  heavily 
doped  GaAs  structures.  Previously,  on  an  RADC  project,  we  had  evaluated 
purchased  GaAs  wafers  that  had  been  doped  during  growth.  However,  the  as 
grown  GaAs  could  not  be  doped  sufficiently  to  achieve  the  high  carrier 
concentrations  necessary  for  the  observation  deformation  potential  coupled 
temperature  effects  and  accordingly  this  study  of  diffusion  in  GaAs  was 
initiated. 

In  order  to  evaluate  the  diffused  layers,  it  was  necessary  to  employ 
selected  area  diffusions  and  surface  acoustic  wave  measurements.  The 
project  then  became  one  of  studying  masked  deep  zinc  diffusions  in  GaAs. 

3.2  Results  and  Conclusions 

The  full  details  of  the  work  are  given  in  the  MSEE  thesis  of  G.  Tuttle. 
The  following  summary  is  taken  from  the  thesis. 

A  capsule  diffusion  system  for  diffusing  zinc  into  GaAs  has  been  built 
and  carefully  characterized.  The  system  satisfactorily  forms  heavily  doped 
zinc  layers  in  GaAs.  Such  deep,  heavily  doped  layers  are  not  possible  with 
ion  implantation.  In  addition,  it  has  been  shown  that  reactively  sputtered 
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silicon  nitride,  in  conjunction  standard  wet  chemistry  processing,  serves  as 
an  adequate  masking  material  to  confine  the  zinc  diffusions  to  selected  area 
of  a  GaAs  wafer.  Replacement  of  the  wet  etching  with  plasma  etching  to 
eliminate  the  anomalous  lateral  diffusion  effects  observed  will  make  silicon 
nitride  a  truly  excellent  masking  system. 

Also,  using  the  capsule  system,  several  attempts  were  made  at 
performing  n-type  diffusions  using  tin  and  sulfur  as  the  donor  material. 
Though  unsuccessful,  these  experiments  provided  some  insight  into  what  type 
of  diffusion  system  might  be  used  for  performing  n-type  diffusion  in  GaAs. 

The  question  naturally  arises  as  to  what  value  the  zinc  diffusions  are 
without  the  n-type  diffusions  needed  for  fabricating  high  performance  GaAs 
transistors.  Perhaps  the  most  obvious  use  is  in  making  p+  contacts  to 
electronic  devices  having  p-layers  as  part  of  their  structure.  In  addition, 
other  semiconductor  devices  such  as  lasers,  LEDs,  and  GaAs  solar  cells  need 
diffused  layers  and  ohmic  contacts. 

If  GaAs  substrates  are  available  with  an  n-type  doping  of  10-16  cm-3  or 
less,  the  zinc  diffusions  and  silicon  nitride  masking  can  immediately  be 
applied  towards  making  a  self-aligned  metal  insulator  semiconductor  field 
effect  transistor  (MISFET) .  All  that  is  needed  is  to  diffuse  the  source  and 
drain  regions  through  a  SiNx  mask,  deposit  and  pattern  the  gate,  source,  and 
drain  metal  concurrently,  and  an  enhancement-mode  p-channel  MISFET  is 
finished  in  two  simple  steps.  The  SiNx  between  the  source  and  drain  serves 
double  duty  as  a  diffusion  mask  and  gate  dielectric.  Albeit,  this  is  a  low 
frequency  device,  but  it  would  still  have  application  as  a  biasing  element 
or  active  load  in  a  high  frequency  digital  or  microwave  circuit. 

The  zinc  diffusions  and  SiNx  masking  could  also  be  used  in  conjunction 
with  ion  implanted  n-type  layers  as  the  diffused  gate  of  a  JFET.  The  JFET 
should  have  certain  advantages  over  the  MESFET.  The  JFET  should  be  a  higher 
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frequency  device,  with  diffusion  the  pinch-off  voltages  should  be  more 
reproducible  from  device  to  device  and  wafer  to  wafer,  and  the  JFET  should 
be  able  to  tolerate  larger  voltage  swings  at  the  gate. 

Also,  using  the  Gaussian  type  diffusions  that  were  demonstrated,  p- 
channel  MESFETs  and  JFETs  (if  used  with  an  n-type  implanted  gate)  should  be 
possible.  However,  more  work  is  required  in  characterizing  these  Gaussian 
diffusions. 

One  last  application  is  related  to  the  acoustic  properties  of  GaAs.  It 
has  been  shown  both  theoretically  and  experimentally  that  the  temperature 
dependence  of  the  elastic  constants  of  silicon  changes  when  the  silicon  is 
heavily  doped  with  boron.  This  can  be  used  to  form  signal  processing 
devices  hat  do  not  change  characteristics  as  the  temperature  changes.  It 
would  be  very  desirable  to  do  this  with  GaAs.  In  ongoing  experiments,  the 
temperature  dependence  of  both  doped  and  undoped  GaAs  is  being  measured 
using  surface  acoustic  waves  traveling  in  undoped  and  heavily  zinc-doped 
GaAs.  If  GaAs  shows  the  same  change  in  temperature  dependence  that  silicon 
does,  temperature  compensated  acoustic  devices  on  GaAs  will  become  possible. 

1.  F.S.  Kahn  and  P.B.  Allen,  "Temperature  Dependence  of  the  Elastic 
Constants  of  p+  Silicon,"  Phys.  Stat.  Sol.(b)  128,  31  (1985). 

2.  J.S.  Wang  and  K.M.  Lakin,  "Low-Temperature  Coefficient  Bulk  Acoustic  Wave 
Composite  Resonators,"  Appl.  Phys.  Lett.  40(4),  15  Feb.  1982. 

3.  G.  Tuttle,  "A  Technique  for  Silicon  Nitride  and  Aluminum  Nitride  Masked 
Deep  Zinc  Diffusions  in  Gallium  Arsenide,"  Iowa  State  University,  1985. 
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4.0  FINITE  DIFFERENCE  CALCULATIONS 

4.1  Introduction 

In  order  to  treat  a  general  class  of  piezoelectric  excitation 
problems,  a  finite  difference  formulation  has  been  developed  to  study  the 
two  dimensional  coupled  piezoelectric  wave  equation  [3].  The  formulation 
allows  for  three  displacement  components  and  an  electric  potential, 
complicated  boundary  configurations,  and  arbitrary  material  anisotropy. 
This  study  will  test  the  formulation  on  a  specific  resonator  structure. 

Some  simple  geometrical  structures  have  been  analyzed  by  others 
using  different  methods  such  as  plate  wave  mode  expansions.  These 
methods  however,  assume  infinite  dimensions  in  some  particular  directions 
and/or  have  restricted  boundary  conditions  [6,7,9,10,  and  12).  They  can 
not  account  for  changes  in  material  thickness  and  require  material 
uniformity  throughout  the  region. 

4.2  Development  of  the  Numerical  Method 

The  theoretical  development  of  the  finite  difference  formalism 
starts  from  the  general  point  form  linear  elastic  equations  for  arbitrary 
anisotropy. 
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du 


k  _  € 


^1  ei,kl  JXj  'il  8xj 


(4.4) 


For  resonator  analysis,  the  time  harmonic  assumption  Is  employed,  as  In 
the  one  dimensional  case,  although  there  are  numerical  methods  for 
solving  time  dependent  wave  equations  as  well. 

The  two  dimension  plane  will  be  the  x^,  X2  plane  with  x^  normal 
to  it.  No  spatial  variations  are  allowed  along  x^.  Equation  (4.1)  can 
then  be  rewritten. 


aXl  Tlj  "  ax,  (_T2j)  “  _P  u  Uj 


(4.5) 


The  left  hand  side  of  Eq.  (4.5)  can  be  written  as  a  cross  product  by 
defining  a  vector  0  where 


0(j)  -  Oj^Xj  +  02(;i)*2 


(4.6) 


and 


0  (J)  •  -t 

1  2j 


°2(J)  "  Tlj  * 


Then, 


Vx  0 


- 


-  pw  Uj  . 


(4.7) 


Equation  (4.7)  can  be  written  in  the  Xj ,  Xj  plane  as 


§  Vx  0^»dS  -  -  u2Jpu  *dS 
•  J 


(4.8) 
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where 


dS  “  dx^dXj 


The  result  of  having  the  left  hand  side  of  Eq.  (4.8)  in  the  fora  of  a  curl 
operation,  allows  a  direct  integration  using  Stoke's  theorem.  This 
reduces  the  expression  by  one  order  to  a  line  integral  of  the  stress 
components  around  a  closed  path  in  the  1,  2  plane, 


$  O^^dl  - 


-  u;  ^pujdx^dx^ 


(4.9) 


Physically,  those  integrals  correspond  to  finding  the  net  force  on  the 
material  region. 

An  expression  for  0^  and  0^^  can  be  found  by  expanding  Eq. 
(4.3)  using  the  corresponding  subscripts, 


n  0>  _  -t  -  -r  ^  r  !!li  .  11  .  11 

1  *2j  U2jkl  dXj  U2jk2  3x2  ejk,l  bx{  ~  e2j,2  3x2 


and 


,<J)  -  -T, 


Su. 


lj  "  Lljkl  aXl  +  Cljk2  3x2  +  elj,l  Sxj  elj ,2  3x2  * 


»“k 


a? 


11 


The  above  equations  can  be  reduced  to 


and 


i  V  ax1  °jk  3x2 


(4.10) 


02<J>  -  -a 


Jk  jk  3x2 


(4.11) 
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by  letting 


aJk " 

C2Jkl 

k-1.2,3 

aj>c“ 

e2J,l 

k-4 

bjk  ■ 

C2jk2 

k-1,2,3 

v  - 

e2j  ,2 

k-4 

Cjk“ 

Cljkl 

k-1,2,3 

Cjk  " 

elj,l 

k-4 

V  “ 

Cljk2 

k-1,2,3 

v- 

®lj,2 

k-4 

1,2,3 

and  u,  - 
4 


Similarly,  Eq.  (4.2)  can  be  rewritten  as 


Di  "  a*2  (~V  "  °  * 


By  defining  a  vector  P  to  be 


P  -  P,*,  ♦  P2*2 


where  Pj  *  -D^  and  P2  “Dp 

Equation  (4.12)  can  be  written  as  the  cross  product 
(  7  x  P)  •  0 

and  using  Stoke* s  theorem  as 
$  P  •  di  •  0  . 


(4.12) 


(4.13) 
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The  expressions  Pj  and  P2  can  he  expanded  with  the  use  of  Eq.  (4.4), 


n  3uk  8  V  duk  3  Y 

°2  "  ”e2  ,kl  3^  21  3^  “  e2  ,k2  ^  +  22  ^ 


(4.14) 


P,  -  D 


du. 


-  € 


a  « 


3u,. 


+  e, 


1  ci,k.i  aXl  11  aXj  1 » w.2  3x2 


-  € 


a  1 

12  ax2  * 


(4.15) 


Reducing  the  above  equations, 


3uk  3uk 

^1  3jk  3x^  ^jk  3x2 


(4.16) 


duk  duk 

P«  *  c ,  *  •.  +  d.,  -5 — 

2  jk  3x^  jk  ax2 


(4.17) 


where 


jk  “  e2,kl 

k-1,2,3 

‘jk  *  f2l 

k-4 

'jk  "  e2,k2 

k-1,2,3 

'jk  "  *22 

k-4 

’jk  "  el,kl 

k-1,2,3 

:jk  “  "*11 

k-4 

‘jk  “  el,k2 

k-1,2,3 

'jk  ”  ~*12 

k-4 

-4“  y* 

The  subscript  arrangement  of  Eqs.  (4.10),  (4.11),  (4.14),  and  (4.15) 
allows  those  four  equations  to  be  written  as 


'1 


•11.  i 


jk  dx,  Jk  Ox, 


(4.18) 
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0<j>  ±k 

C.i  >  '  o  . .  \ 


'Jk  dx 


jk  3x. 


(A. 19) 


where  j-1,2,3,4  and  k-1,2,3,4, 
and  che  cocal  vector  Q  Is 


q  Cj)  „  +  Q2(;i)x2  , 


which  represents  the  sura  of  Eqs.  (4.10),  (4.11),  (4.14),  and  (4.15).  The 
material  censors  ajjt»  h^*  cjk*  anc*  ^jk  are  aS  Previous^y  defined. 

Thus,  there  are  four  line  integrals  of  the  form 


£  Q  •  dl  -  J’Q?dx2  +  /o^dx^  ,  (4.20) 

composed  of  16  terms  in  each  of  4  regions  for  a  total  of  256  integrals, 

and  three  area  integrals  from  Eq.  (4.9)  since  the  right  hand  side  of  Eq. 

(4.13)  is  zero. 

To  evaluate  the  line  integral  given  by  Eq.,(4.20)  consider  the 
geometry  shown  in  Fig.  4.1.  The  right  hand  side  of  the  equation  can  be 
written  as 

b  c  d  e  f 

•  dl  ■  /  02dx2  +  S  QjdXj  +  J  QjdXj  +  J  Q2dx2  +  S  Q2&X2 

abode 


5  h  a 

+  J  Q.dx.  +  /  QjdXj  +  /  Q2dx2 

f  g  h 


(4.21) 


where  the  subscripts  have  been  dropped  and  the  integration  limits  refer 
to  the  various  line  segments  in  che  four  material  regions,  denoted  by  0  , 
(2)  »  Q)  *  and  © ,  of  Fig.  4.1. 


mm 


Fig.  4.1.  A  rectangular  sampling  grid  in  the  1,2  plane  showing  the 

method  of  labeling  mesh  points  and  material  regions.  The  line 
Integral  is  carried  out  along  the  a-b-c-d-e-f-g-h-a  path  and 
the  area  integral  within  Che  region  defined  by  the  line.  The 
encircled  numbers  refer  to  a  particular  quadrant  and  its 
corresponding  material  constants  (p,«,  C,  e) 
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Using  Eqs.  (4.18)  and  (4.19)  che  right  hand  side  of  Eq.  (4.21)  can 
be  evaluated  and  written  as  a  function  of  a,  b,  c,  and  d,  and  the 
particle  displacement  and  potential  fields.  For  notation  simplification 
purposes,  the  mesh  points  will  be  referred  to  as: 

(I.J)  0 

(l.J+l)  l 
(I+l.J)  2 

(I.J-l)  3 

(I-I.j)  4  . 


The  terms  involving  Qn  will  be  evaluated  first, 


S  Q2<!x2  -  S  Pf;  <1*2  *  S  ^  d*2 

a  a  a 


P(  }  ♦  PU(M  . 

nl  L  0 


-  (u(l )  -  u(0))  +  P(  u(2)  -"--U-°--  )  (4.22) 


where,  in  general. 


du  m  u(l)  -  u(0) 
dx  h. 


i  -  1-4, 


and  because  the  points  are  sufficiently  closed  together. 


u(b)  -  piLLl 


mm 


mmm 


In  a  similar 

way, 

a 

/  Q,dx 
h-  z 

.  c©(  - 
h 

u(0 )  . 

)  x, 

1  L 

|  +  <^u(a: 

-h4/2 

&  ,  u(0)  -  u(4 

2hi 

(  2 

u(h)  - 

u(4)  +  u(l) 

2 

• 

X  Q2d,3 

,  -  c®(  u(0> 

-  u(3)  x 

)  x 

,  |  +  &  (u(e 

h2/2 

Ph 

-  (u(0)  -  u(3))  -Jj-i 

C-l)  +  P(  u(0) 

u(e)  ■ 

u(3)  +  u(0) 

2 

u(3)  +  u(2) 

(4.27) 


)  -  t^(u(2)  -  u(0 ) )  (yj±) 


where 


u(g) 


u(4)  +  u(0) 


S  Q1dx1  -  -P  (u(g)  -  u( f ) )  -  tP(  ..:..u.(4).  )  Xl  I 


-h3/2 


_aG>  (  h<o).,.- ' all!  )  -  P  (u(o)  -  u(4 ) )  . 

2  2h4 


And  finally. 


J  QjdXj 


-£>c„(h)  -  u(g);  -  g>(  ,  %1  |”‘ 


hx/2 

0 


-P  (  ~  U-°--  )  -  b®  (u(0)  -  u(4))  yy  . 

4 


The  total  line  integral  can  be  found  by  adding  Eqs.  (4.22) 
(4.29).  The  sum  can  be  simplified  to 


•  dl  -  (u(l)  -  u(0))  +  (u(2)  -  u(0))a2 


(4.28: 


(4.29: 


through 


+  (u( 3 )  -  u(0) )  a  +  (u(4)  -  u(0) )  a,. 


(4.30 
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A,  *  ft,  m 

'i  *  ,<  ‘ 2h^  3  >  *  )  1  . 


ft4  ♦  ft2  ,®_  P 
V  2h3  +  <^-^>1  . 


ft3  *  bft  P_  P 

34-l‘— 5h7 - , 

4 


where  the  subscript  In  1  refers  to  a  particular  mesh  point. 


Consider  now  the  area  integral  suggested  by  the  right  hand  side  of  Eq. 


(4.9), 


-■£  J'ou(0)dx1dx2  »  -  <j?  u(0)J*  pdx^x2  . 


The  evaluation  of  the  are  integral  yields 


where 


-o?  u(0 )j  pdx1dx2  -  -  Q0u(0) 


Qq  -  of?  , 


(4.31) 


where 


p  ■  lp^i1h2  +  f^2h3  +  +  P^4hi^ 


is  the  effective  mass  per  unit  depth. 

Since  Eqs.  (4.30)  and  (4.31)  are  an  expanded  form  of  the  left  and 
right  hand  side  of  Eqs.  (4.9)  and  (4.13),  they  can  be  set  equal  to  each 
other.  After  the  subscripts  are  restored,  the  expression  becomes 
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(\<‘>  *  “kw»*ijk  ♦  <“k<2>  -  ♦  <\‘3>  -  \<°»V 

♦  (uk(«)  -  <^(0))!^  ■  Q0u,(0). 


where 


k,j  -  1.2. 3, 4 


QQ  «  0  for  j  ■  4, 


che  repeaCed  subscript  indicate  a  sum,  and  the  first  subscript  in  Cl 
refers  to  a  particular  mesh  point. 

The  above  expression  can  be  rewritten  as 


where 


uk(°)ajk  -  QqUj (0)  -  Uj 


Uj  “  uV.(l)aljk  +  “k^^jk  +  “k{3)a3jk  +  “k^^jk 


ajk  ^Ijk  +  a2jk  +  a3jk  +  a4jk  • 


(A. 32) 


and 


P  h  ♦  c®  h  P  -  .© 

jk  2  Jk  4  +  Jk  a  jk 


ljk  2h 


a  _  ^Jkhl  *  ^jkh3  ^jk  '  ^Jk 
2jk  2h,  2 


P  h  +  c®  h  a© 

Jkh4 _  Jkh2  jjk. 


3jk  2h 


JJi 
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a  .  jjkh  *  fjkh  ,  Pj>  - 

4jk  2h4  2 


Equation  (4,32)  holds  at  each  point  in  the  mesh  and  can  be  written 
as  a  matrix  as  follows, 


p 

K- 

1 

o 

o 

a!2 

*13 

1 

1 

1 

1 

*14 

r^w) 

V 

*21 

<*22-0 

o 

*23 

1 

1 

1 

*24 

j u2(0) 

U2 

*31 

a32 

a33'Qo 

1 

1 

34 

ju3(0) 

U3 

a41 

m 

a42 

*43 

1 

1 

1 

a44_ 

[u4(0) 

L°‘J 

(4.33) 


where  u(0)^_3  stand  for  the  particle  displacement  fields  corresponding  to 
the  central  mesh  point  along  the  direction  given  by  the  subscripts,  and 
refers  to  the  electric  potential  corresponding  to  that  same  central 
mesh  point.  Physically,  the  fields  at  a  given  mesh  point  are  viewed  as 
being  driven  by  the  nearest  neighbors.  The  dashed  lines  in  the  matrix 
separate  those  parts  coupled  by  piezoelectricity. 

The  material  tensors  a,b,c,d,  can  be  written  in  matrix  form  as 
follows. 


V 

k-1 

2 

3 

4 

J-l 

C16 

C66 

C56  : 

1 

*16 

2 

C12 

C26 

C25  j 

*12 

3 

C14 

C46 

C45  ; 

14 

e 

4 

e12 

*26 

e25  ! 

-  f21 

(4.34) 


J-l 


1  *26 

k-1 

*22 

2 

*24 

3 

C11 

C16 

C15 

C16 

C15 

C66 

C56 

C56 

C55 

*11 

*16 

*15 

I 

2 

3 

Each  censor  is  given  in  terms  of  the  usual  reduced  subscript  material 
tensors  and  are  partitioned  by  the  dotted  lines  to  identify  those  regions 
associated  with  Laplace's  equation  only,  the  wave  equation  without 
piezoelectric  coupling,  and  finally  the  piezoelectric  coupling 
coefficients  [3]. 

With  the  proper  simplifications  the  two  dimensional  method  can  be 
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used  to  do  one  dimensional  analysis. 

4.3  Application  of  the  2D  Method 

The  two  dimension  numerical  method  will  be  applied  to  a  structure 
consisting  of  an  electrode,  a  piezoelectric  thin  film  layer,  and  a 
semiconductor  substrate  (see  Fig.  4.2).  The  mass  and  the  thickness  of 
the  electrode  are  going  to  be  considered  negligible  in  this  test  of  the 
formulation.  The  energy  trapping  effect  will  then  be  solely  due  to  the 
finite  dimensions  of  the  structure.  A  standing  wave  pattern  due  to  the 
material  discontinuities  along  the  two  dimensions  under  consideration 
will  be  created  within  the  acoustic  cavity  when  an  electric  potential 
field  is  applied  to  it. 

The  particle  displacement  field  perpendicular  to  the  1,2  plane  will 
be  considered  uncoupled  from  the  electric  potential  field  by  assumption. 
This  Implies  that  the  (3,4)  and  (4,3)  positions  in  the  tensor  of  Eqs. 
(4.34)-(4.37)  will  be  set  to  zero.  Three  fields  will  be  considered  in 
the  piezoelectric  region:  a  potential  field  and  two  particle  displacement 
fields,  one  along  each  dimension.  In  the  substrate  region,  only  the  two 
particle  displacement  fields  will  be  considered  since  the  potential  field 
decouples  because  of  the  material  properties  of  the  substrates  used. 

For  simplicity,  the  electric  fields  fringing  out  the  resonator 
structure  will  not  be  included  in  this  work.  It  was  considered  that  with 
an  8:1  dielectric  ratio  in  favor  of  the  piezoelectric  material,  the 
electric  fields  escaping  to  the  outside  of  the  cavity  region  could  be 
neglected.  To  achieve  the  nonfringing  effect  of  the  stress  field,  the 
material  constant  p,  c,  e  will  be  set  to  zero.  This  forces  the 
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continuous  across  the  boundary.  Fringing  could  be  accounted  for  by 
considering  a  mesh  larger  than  the  acoustic  cavity  structure.  Jor 
continuity  reasons  also,  the  electric  potential  field  and  the  force  oust 
be  equal  along  the  '2*  direction  across  the  nonpiezoelectric- 
piezoelectric  boundary.  The  value  of  the  electric  potential  at  the 
electrodes,  a  known  independent  variable,  becomes  a  part  of  Eq.  (A. 33) 

(see  Fig.  A. 3). 

The  resonator  structure  of  interest  requires  chat  the  piezoelectric 
crystal  C-axis  be  oriented  along  the  '2'  direction.  The  material  tensor 
matrices  for  the  piezoelectric  material  will  then  have  to  be  transformed 
to  match  the  coordinate  axis  in  use.  The  transformation,  see  Fig.  A. A, 
is  simpler  to  realize  when  the  material  tensor  elements  are  written  using 
the  abbreviated  subscript  notation. 

As  previously  mentioned,  the  resonator  structure  will  be  divided  by 
a  mesh  of  vertical  and  horizontal  lines.  The  value  of  the  electric 
potential  and  particle  displacement  field  for  each  mesh  point  will  be 
found  by  using  a  modified  version  of  Eq.  (A. 33).  The  third  column  and  row 
of  the  system  in  Eq.  (A. 33)  will  be  deleted  since  the  particle  displacement 
field  along  the  '3'  direction  is  considered  uncoupled  as  mentioned.  When 
expanded  through  the  mesh  structure,  the  sec  of  difference  equations  can 
be  written  in  the  form  of  a  sparse  banded  matrix  similar  to  the  one 
obtained  for  the  one  dimensional  case. 

For  programming  purposes,  the  frequency  dependent  expression  Q  will 

o 

be  normalized  to  the  fundamental  parallel  frequency  fQ  of  the  thickness 
longitudinal  mode  in  an  all  Si  structure.  The  expression  for  Qq  was 


numerical  results  will  be  qualitatively  evaluated. 

As  was  noted  in  chapter  III,  a  certain  minimum  number  of  mesh  points 
is  needed  to  determine  the  field  waves  with  accuracy.  Extremely  large 
matrices  can  be  generated  if  the  same  resolution  is  desired  for  the  two 
dimensional  case.  For  this  reason,  the  modeling  will  be  limited  to 
frequencies  close  to  fundamental  where  the  fields  can  still  be  resolved 
using  a  small  number  of  mesh  points.. 

The  particle  displacement  at  a  particular  mesh  point  will  be 
represented  by  a  dot.  For  clarity  purposes,  only  the  waves  along  the 
perimeter  of  the  deformed  structure  will  be  shown  with  solid  lines.  The 
undeformed  mesh  and  device  geometries  are  included  as  reference.  The 
electric  potential  field  will  only  be  shown  in  the  piezoelectric  region 
since  it  is  zero  in  the  rest  of  the  structure. 

The  particle  displacement  magnitudes  shown  in  Figs.  4. 5-4. 7,  are  by 
no  means  proportional  to  the  dimensions  of  the  structure.  Even  though 
the  latter  are  normalized,  they  are  generally  in  the  micron  unit  range 
whereas  the  particle  disp1 aceraents  are  in  the  Angstrom  unit  range. 

Physical  considerations  suggest  that,  due  to  the  structure  of  the 
crystals  used,  the  particle  displacement  field  along  the  '1'  direction 
and  the  magnitude  of  the  electric  potential  field  be  symmetric  with 
respect  to  a  vertical  axis  along  the  middle  of  the  acoustic  cavity.  It 
is  expected  that,  because  of  the  Poisson's  ratio^,  the  particle 

1 1 1  refers  to  the  ratio  of  the  transverse  contraction  per  unit 
dimension  of  a  bar  to  its  elongation  per  unit  length,  when  subjected  to  a 
tensile  stress. 


displacement  along  the  *  1 '  direction  be  larger  in  the  ZnO.  If  the 
particle  displacement  and  potential  field  results  are  correct,  V»D  -  0. 
Also,  u  and  ¥  should  be  continuous  across  the  ZnO/Sl  interface.* 

The  expected  field  symmetry  can  be  seen  in  Figs.  4. 5-4. 7.  The 
narrow  structure  of  Fig.  4.5a)  clearly  shows  the  effect  of  the  Poisson 
terms  as  the  structure  deforms  along  a  direction  different  than  that  of 
the  driving  electric  potential.  It  can  also  be  seen  from  this  figure 
that  the  particle  displacement  along  the  *  1 '  direction  is  larger  in  the 
ZnO  region.  The  structure  however,  is  not  wide  enough  to  sustain  a 
standing  wave  resonance.  Such  a  wave  pattern  can  be  observed  in  Figs. 
4.6a)  and  4.7a)  as  the  horizontal  to  vertical  dimensions  ratio  of  the 
structures  gets  larger.  Inversely  proportional  to  this  ratio  is  the 
effect  of  the  Poisson  terms.  The  electric  potential  is  also  symmetric  as 
shown  in  Figs.  4.5b)-4.7b).  The  electric  potential  variations  along  the 
x^  direction  are  directly  related  to  the  magnitude  changes  in  the 
particle  displacement  field.  These  variations  however,  are  small  and 
impossible  to  resolve  in  the  figures. 

A  remark  about  the  behavior  of  the  particle  displacement  field  in 
the  Fig.  4.5  structure:  the  larger  particle  displacements  in  the 
piezoelectric  region  are  not  only  due  to  the  Poisson  ratio  but  also  to  an 
acoustic  impedance  mismatch.  It  happens  that  a  large  part  of  the  wave 
generated  in  the  piezoelectric  region  is  reflected  at  the  ZnO/Si 
boundary.  The  mismatch  decreases  as  Che  structure  becomes  wider  as  can 
be  seen  from  the  behavior  of  the  field  in  Figs.  4.6  and  4.7. 

The  no  free  charge  condition  was  tested  by  applying  the  equation 
•D  -  0  to  a  mesh  point  in  the  piezoelectric  region.  As  expected,  the 


Fig.  4.6.  Field  configuration  In  a  12  vertical  by  17  horizontal  me ah  point  structure  pt;  an 

effective  normalized  frequency  of  1.27.  The  horizontal  to  vertical  spacing  ratio  le 
a)  particle  displacement  fluid  ( 1  cm  -  1  A  displacement ,  b)  electric  potential  field 
(volts),  where  the  vertical  lines  represent  0  V  reference  lines 
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Fig.  4.7.  Field  configuration  in  a  20  vertical  by  15  horizontal  mesh 

point  structure  at  an  effective  normalized  frequency  of  1.27, 
The  horizontal  to  vertical  spacing  ratio  is  2.  a) 
displacement  field  ( 1  cm  ■  1  X  displacement),  electric 
potential  field  (volts) 


condition  holds.  The  continuity  of  the  particle  displacement  field 
across  the  boundary  can  be  visually  confirmed  by  looking  at  the  figures. 
The  waves  along  both  directions  seem  to  be  continuous  across  the  material 
interface.  A  qualitative  confirmation  of  particle  displacement 
continuity  is  affected  by  the  number  of  mesh  points  required.  The 
continuity  of  the  electric  potential  field  across  the  discontinuity  is 
obvious.  The  electric  potential  is  zero  at  the  interface,  in  this 
analysis,  and  so  is  the  electric  potential  in  the  nonpiezoelectric 
region. 

An  impedance  expression  can  be  found  by  considering  the  electrical 
properties  of  che  structure  at  the  terminals.  As  for  the  one  dimensional 
case,  consider  the  expression 


Z 


V 

iujDa 


(4.41) 


For  the  two  dimensional  case,  the  electric  displacement  D  at  the 
electrode  can  be  approximated  by 

L-l 

j’D  •  da  »  a  Z  D„(n) 
n*l 


where  a  is  the  electrode  area  and  L  is  the  number  of  mesh  points  along 
the  *1'  direction.  The  electric  displacement  D^,  along  the  '2' 
direction,  is  given  by 
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where  the  material  tensor  transformation  has  already  been  performed.  The 
above  equation  reduces  to 


L-l  L-l  Buj  du^ 

“  °2  “  (a33  1Z  ~  *33  TT  )m 

n-1  n-1  2  2 


(4.43) 


The  impedance  expression  can  be  found  by  substituting  Eq.  (4.43)  into 


Eq.  (4.41), 


L-l  e.«  Su«  du, 

jwC  d_  (  £  (  -  — —  )) 

g  2  n-1  *33  3x2  dx2 


(4.44) 


where 


C«'V 


i  •  -i  . 


a  is  the  area  of  the  electrode  in  the  Xj,  plants,  and  d 2  is  the 
thickness  of  the  piezoelectric  region  where  the  capacitance  effect  is 
concentrated. 

For  computational  purposes  the  Impedance  will  be  normalized  to  JwC  . 

S 

Figure  4.8  shows  a  plot  of  the  magnitude  and  phase  of  the  impedance  at 
frequencies  around  fundamental  series  and  parallel  resonances.  Only  one 
resonance  spectrum  can  be  observed  since  just  one  mode  is  trapped  in  the 
structure  of  12  vertical  and  5  horizontal  mesh  points,  a  structure 
similar  to  the  one  shown  in  Fig.  4.5.  Computational  considerations  were 
again  a  factor  in  choosing  those  particular  dimensions. 


Fig.  4.8.  Magnitude  and  phase  Impedance  plot  of  a  12  vertical  by  5 
horizontal  mesh  points  structure  at  frequencies  around 
fundamental  series  and  parallel  resonance.  The  materials 
intersect  at  the  sixth  vertical  mesh  point  from  the  bottom 


At  frequencies  below  resonance,  the  normalized  impedance  tends  to  1 
as  in  the  one  dimensional  case.  Note  that  the  fundamental  parallel 
resonance  will  not  occur  at  a  normalized  frequency  of  1,  becautfe  of  the 
previous  normalization  of  the  frequency  to  a  parallel  resonance  frequency 
of  an  all  Si  acoustic  cavity.  The  overtones  though,  will  still  occur  at 
approximate  integer  multiplies  of  the  fundamental.  The  series  resonance 
will  vary  as  given  by  Eq.  3.24. 

The  numerical  impedance  results  and  Eq.  3.24  can  be  used  to 
calculate  an  effective  coupling  coefficient  K4, 


can 

Ts 


2 

from  where  K  ■  0.062.  This  value  can  be  compared  to  one  calculated  from 
the  elastic,  dielectric  and  piezoelectric  coefficients  of  ZnO, 


where  m  * 


'33 


CE  €S 
C33  33 


and  m 


15 


CE  eS 

C55  11 


,  for  a  longitudinal  wave, 


for  a  shear  wave. 


The  effective  coupling  calculated  from  the  impedance  is  35Z  smaller  than 
the  analytic  coupling  calculated  for  a  longitudinal  wave,  and  8%  larger 
than  the  analytic  coupling  calculated  for  a  shear  wave.  And  therefore  it 
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Is  concluded  Chat  Che  effective  coupling  coefficient  calculated  from  the 
Impedance  Is  a  combination  of  the  coupling  coefficients  calculated  from 
the  two  waves. 


4. 
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4.4  Conclusion 

A  two  dimensional  difference  formulation  of  the  general  anisotropic 
piezoelectric  coupled  wave  equations  has  been  implemented  for  the  modeling 
of  the  bulk  acoustic  wave  resonator's  wave  structure.  The  difference 
equations  form  a  coupled  set  of  four  equations  describing,  in  a  straight 
forward  manner,  the  particle  displacement  and  potential  fields  at  each 
element  in  the  mesh.  The  form  of  the  equations  indicate  that  a  central 
mesh  point  is  driven  by  its  nearest  neighbors.  When  expanded,  the 
difference  equations  can  be  written  in  the  form  of  a  sparse  banded  matrix 
whose  solution  is  obtainable  by  well-known  techniques. 

A  one  dimensional  comparison  of  analytical  and  numerical  solutions 
indicates  that  the  finite  difference  method  can  be  applied  with  excellent 
results. 

The  two  dimensional  numerical  results  satisfy  the  requirements  of  no 
free  charge  and  symmetry  with  respect  to  a  vertical  line  across  the 
middle  of  the  structure.  The  continuity  of  the  particle  displacement  and 
potential  field  across  the  ZnO/Si  interface  is  visually  confirmed.  A 
quantitative  check  was  not  carried  out  because  the  number  of  mesh  points 
required  to  do  so.  A  standing  wave  resonance  along  the  x^  direction  is 
observed  when  the  horizontal  dimension  is  larger  than  half  a  wavelength. 

A  series  and  parallel  resonance  confirm  that  the  structure  can  resonate. 

A  comparison  with  analytical  results  obtained  by  others,  was  not  possible 
because  of  the  difference  in  analyzed  structures. 

The  two  and  one  dimensional  analysis  results  suggest  that  the 


46 


numerical  method  can  be  successfully  used  in  modeling  the  two  dimensional 
pattern  inside  a  bulk  acoustic  wave  composite  resonator. 


4  W 
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